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The two-band Hubbard model with the density of states obtained from the band calculation 
is applied for FeSi, which is suggested to be a Kondo insulator or a correlated band insulator. 
Using this model, the correlation effects on FeSi are investigated in terms of the self-consistent 
second-order perturbation theory combined with the local approximation. The calculated opti- 
cal conductivity spectrum reproduces the experiments by Damascelli et al. semiquantitatively 
and the specific heat explains the anomalous contribution at about 250 K observed in FeSi. 
Inclusion of the spin fluctuation and the extension to the case of strong correlation are also 
discussed. 
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In the study of a specific material among the strongly correlated electron systems, the effect 
of the band structures often plays a crucial role when one compares a theoretical calculation to 
the experiments. Use of a simple theoretical model might not capture the salient features of the 
material. Development of a theoretical method that is capable of taking proper account of the 
realistic features of the material is necessary. We report our recent approach to the study of the 
anomalous properties of FeSi in such direction. 

FeSi is well known for more than thirty years and a number of studies from various aspects have 
been done, stimulated by the fascinating physical properties. The early study by Jaccarino et al% 
showed that the susceptibility is much enhanced over the value expected from the band param- 
agnetism at finite temperatures and has a broad peak at about 500 K. It was also reported that 
the specific heat seems to have an anomalous enhancement at about 250 K. These behaviors were 
explained by a band model with an energy gap, but unphysically narrow bands were necessary, so 
that this difficulty has attracted interests of many researchers. From the conductivity measure- 
ments, FeSi is an insulator at low temperatures but shows metallic behavior at room temperature. 
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** e-mail address: saso@phy.saitama-u.ac.jp 
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To explain these unusual properties of FeSi, several theoretical approaches have been proposed, but 
the most successful one is the spin fluctuation scenario by Takahashi and Moriya.& It explains the 
anomalous magnetic property of FeSi and their idea of the thermally induced magnetic moment 
was confirmed by the neutron scattering experiment 

The recent optical studies,Blli'& however, revealed the unusual properties of FeSi again. 
Schlesinger et al. reported that the gap of about 60 meV (~700 K) opened at low tempera- 
tures is filled and almost closed at room temperature (about 250~300 K), which they attributed 
to the correlation effect. The following experiments also reported the evidence of the correlation 
effects at low temperatures.i'i'0'0^ In these contexts, Aeppli and Fisk0 suggested that FeSi can 
be viewed as a Kondo insulator or a strongly correlated insulator. 

Kondo insulators have been found in the f-electron systems and typical examples are YbB 12 © and 
Ce 3 Bi 4 Pt 3 © and so on. They have correlated f-bands and small energy gaps at low temperatures. 
Although there are many similarities among FeSi and these materials, the correlation in FeSi 
may not be so strong. However, the same physics can be recognized both in FeSi and Kondo 
insulators, if one reexamines the experimental data carefully. From this aspect, Fu and Doniach© 
proposed an extended Hubbard model with two mixed conduction bands, which is based on their 
band calculation© for FeSi, and confirmed the importance of the correlation effects in physical 
quantities. Their calculation, however, seems to include some errors about the treatment of the 
self-energies. Therefore, we reinvestigated this model carefully and calculated the correlation effects 
in more correct way© and confirmed that the correlation effects do play important roles, but the 
shape of the spectrum in the optical conductivity did not coincide with the experimental data, 
because of the use of the too simple model Hamiltonian. 

Therefore in the present report, we use an extended two-band Hubbard model with the density of 
states obtained from the band calculation, and attempt to explain the low temperature anomalies 
of FeSi observed in the optical conductivity© and the specific heat consistently. 

The band calculations© QBIEil' for FeSi predict that the ground state is a band insulator and 
a recent calculation© reproduces the gap size close to the observed one. Therefore, we start from 
the band insulator model, which consists of two Hubbard bands for d-electrons as follows. 
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where the c| (cj ao -) creates (destroys) an electron on site % in band a =1, 2 with spin a. The tight 
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binding parameters i?- should be fitted to the band calculation and U, U2, U3 and J denote the 
Coulomb and exchange interactions. 

Since one can expect that the optical conductivity spectrum reflects the structure of the quasi- 
particle density of states (DOS) of a system, we use the DOS obtained from the band calculation 
for FeSi by Yamada et al^ for the initial DOS so as to enable detailed comparison with the 
experiment. 

Furthermore, we start from the following general expression of the current operator, 

j = ^E"k m '4vk, (2) 

o-,k mm' 

where m denotes the band indices and derive the convenient expression for the optical conductivity. 
For simplicity, we set the intra- and interband contributions to be equal [v^ m ' = v^). Moreover, we 
assume that the momentum conservation is violated in real systems by some defects and phonon- 
assisted transitions. Therefore, using the linear response theory, we consider the current-current 
correlation function as below, 
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where T^ff' a (iv\ \ iu> n ) denotes the vertex function, and set [. . .] constant. For the present case, this 
leads to the following expression for the optical conductivity, 

n v J- oa w 

x KH + [Pi (v + u) + p° 2 {y + u)], (4) 

where p^{v) denotes the DOS for the band a. This joint-DOS-like form for the optical conductivity 
is simple but convenient for the present case. We set (ev) 2 /% = 1 for simplicity. 

Firstly, we show the optical conductivity obtained from the Hartree-Fock approximation (HFA) 
or a rigid band model in Fig. 1. The used DOS is displayed in Fig. 2 by the solid line for T = 0. 
The DOS is independent of the temperature within HFA. At K, only the interband contribution 
survives and reproduces the shape of the spectrum of the experiment at 4 K in Fig. 3. Therefore, 
the band calculation by Yamada et al. seems to give a good result about the whole structure of 
the DOS at T = but with a slightly smaller gap size (see the comparison with the experiment 
below). Within the rigid band model, however, since the gap is filled only with the intraband 
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Fig. 1. Fig. 1 Calculated optical conductivity within the Hartree-Fock approximation or a rigid band model. 



(Drude) contribution, the temperature variation is monotonous and the spectrum does not become 
flat at a temperature of the order of the gap size. This disagreement was shown by Fu et al. 
first. Ohta et a3 also calculated the optical conductivity in the joint-DOS form from their band 
calculation, but the flat part of the optical conductivity spectrum within the gap could not be 
reproduced. Therefore, the rigid band model is not sufficient to explain the experiments. 

Next, we investigate the correlation effect in the low energy and low temperature region of this 
model. Therefore we calculate the correlation effect by the self-consistent second-order perturbation 
theory (SCSOPT) combined with the local approximation. The second-order self-energies are given 
by 

sS 2)<T (w) = JJJ°° dsid£ 2 d£ 3 
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where 
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Here, N is the number of sites, /(e) the Fermi function and p^{e) the DOS of band a for the 
non-interacting case. To make numerical calculation easy, we take 5 finite (5 = 1CP 7 ) in eq. (6) 
and convert these equations with the transformations© 



These equations have to be solved self-consistently. In this paper, we set U2 = U—J and U3 = U—2J 
in order to reduce the number of parameters. In this case, the Hamiltonian is rotationally invariant 



In the following results, U = 0.5 eV and J = 0.35C7 are chosen so as to reproduce the shape and 
the temperature dependence of the optical conductivity spectrum. The solid line for T = in Fig. 
2 indicates the initial DOS at K, and the correlation effect is absent except the Hartree-Fock 
contribution since the band 1 is filled and the band 2 is empty. 

Note that the gap in the DOS is widened by 16 % so as to reproduce the shape of the spectrum 
of the optical conductivity at 4 K in the experiment, which does not change the essence of the 
following result. Then, the gap size (E g ) of 75 meV is obtained if the steepest parts of the DOS at 
the both sides of the gap are extrapolated and the tails are neglected. (If we regard the gap as the 
region inside the tails of the gap edge, we obtain 60 meV.) The band 1 and 2 in our Hamiltonian 
correspond to the upper and lower part of the DOS with respect to the Fermi level {Ep = 0) 
as is seen in Fig. 2, where we introduce a cut off for each band so as to include one state per 
spin in each band. Then the band width for the band 1 and 2 are about 0.56 eV and about 0.85 
eV, respectively. Although the DOS is asymmetric, the chemical potential is fixed at u = and 
assumed to be temperature independent. One can see in Fig. 2 that the correlation is introduced 
at finite T through the thermally excited electrons and holes and the gap existing at K is almost 
filled up at the temperature of the order of its size, which results in the temperature variation of 
the interband contribution of the optical conductivity (see below) . 

In Fig. 3(a), the temperature variation of the optical conductivity calculated from the 
temperature-dependent DOS in Fig. 2 is shown. In our calculation (Fig. 3(a)), the gap is al- 
most filled up at 300 K as well as the rapid increase in the gap region from 150 to 300 K is seen. 
This is consistent with the experiment (Fig. 3(b)), where the gap is filled rapidly from 100 K to 
300 K. Reflecting the correlation effects, the peak at the gap edge shifts to lower frequency region, 
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if the two bands are degenerate.Ej) 
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Fig. 2. The temperature dependence of the quasi-particle DOS and the initial DOS obtained from the band cal- 
culation(Ref. 26) at T = 0. At finite T, the DOS is strongly temperature dependent due to the correlation 
effects. 



as is seen in the experiment. In our calculation, however, there are dips between the Drude and the 
interband contributions in contrast to the experiment. This may be caused by the simplification 
in deriving eq. (4). However, the almost flat spectrum is obtained at 300 K, which comes from the 
temperature dependence of the interband contribution. 

We also calculate the temperature variation of the specific heat with the same parameters as in 
the optical conductivity. Starting from the equation of motion,© we obtain the following expression 
for the total energy per site: 



1 r°° 

E = oE / ^/MNtfM + ^M} 

^ „ J -co 



+4E { £ IpKKu) +4P2(k,w)}] > ( 8 ) 

k 

where (e^) is the Fourier transformation of tjj (tfj)- The specific heat can be calculated from the 
numerical differentiation of the energy as Cy = (dE/dT)y . The difference between the cases with 
U = and 0.5 eV in Fig. 4 indicates the contribution from the correlation effect, which results in 
a peak of about 4 J/K mol at about 250 K, and explains the "anomalous" contribution (~6 J/K 
mol) in the specific heat at about 250 K reported by Jaccarino et a/.& Note that they evaluated the 
anomaly by subtracting the specific heat of CoSi after the normal electronic contributions 7FeSi an d 
7CoSi are removed, respectively. In the above calculations, we confirmed that the correlation effect 
is essential to explain the temperature dependence of the optical conductivity and the specific heat 
in FeSi. At higher temperatures or for magnetic properties, however, it is also important to take 
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Fig. 3. (a) The temperature dependence of the optical conductivity calculated with the eq. (4). (b)The experimental 
data from Ref. 11. The peaks due to phonons observed in the gap are omitted. 
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Fig. 4. Calculated specific heat using the same parameter as in Fig. 3(a). 
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the spin fluctuationsBEZP into account. 

The self-consistent renormalization (SCR) theory of spin fluctuations has succeeded in describing 
the itinerant magnetism and the quantum critical phenomena with a small number of parameters.© 
On the other hand, the dynamical mean field theory (DMFT) is one of the most powerful schemes 
to take account of the strong local correlation. One of the authors has proposed a new and practical 
scheme that unifies DMFT and SCR in a microscopic way© Application of this theory to FeSi 
may improve the present calculation towards the inclusion of the effects of spin fluctuations at finite 
temperatures and the intermediate coupling. 
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